Decoherence of Bose-Einstein condensates in microtraps 
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We discuss the impact of thermally excited near fields on the coherent expansion of a condensate in a minia- 
turized electromagnetic trap. Monte Carlo simulations are compared with a kinetic two-component theory and 
indicate that atom interactions can slow down decoherence. This is explained by a simple theory in terms of the 
condensate dynamic structure factor. 
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The decoherence of atomic de Broglie waves is a key issue 
for applications in atom interferometry and quantum informa- 
tion processing. It is particularly relevant for integrated atom 
optics based on miniaturized hybrid electromagnetic surface 
traps 1 1, 2] because the atoms couple to a macroscopic, 'hot' 
substrate nearby. Loss processes due to spin flips driven by 
thermal magnetic near fields have very recently been observed 
in the laboratory |3], in agreement with predictions made by 
one of us |4]. In this paper, we discuss a simple decoherence 
scenario for Bose-Einstein condensed atomic matter waves in 
a quasi-one-dimensional microtrap. This setup provides a re- 
alization of the standard model of environment-induced deco- 
herence |5] featuring two attractive advantages: (i) the cou- 
pling to the environment can be microscopically modelled in 
terms of the magnetic dipole interaction; (ii) due to atomic in- 
teractions, the matter wave equation becomes nonlinear and 
novel features are expected. We compare Monte Carlo simu- 
lations for the condensate order parameter to a kinetic theory 
for the matter wave coherence function and show that already 
for moderate interaction parameters, a Bose-Einstein conden- 
sate is more robust with respect to a fluctuating environment. 

We consider an elongated trap similar to those formed 
above current carrying wires In the confinement dom- 
inated regime, the matter waves can be described in a one- 
dimensional mean field approximation |6] (units with Ti = 
m = 1), 



id t ifj 



v(mW + sMm)IV, id 



where the interaction parameter g = 2£l r a/(l — 1.46a/a r ) 
depends on the three-dimensional scattering length a, the ra- 
dial confinement frequency il r and ground state size a r 
The density |-0(x,i)| 2 is normalized to the total number of 
particles N. The potential V(x, t) determines the dynam- 
ics in the axial direction. We assume that for t < 0, the 
atoms are confined in a harmonic trap with frequency O and 
occupy all the zero-temperature condensate mode 4>q(x) @]. 
For t > 0, the axial confinement is switched off, the atoms 
expand, and we take into account their interaction with mag- 
netic field fluctuations by letting V(x, t) be a random poten- 
tial. Note that the radial confinement is kept constant. Eq. ([0 
thus describes the interplay between matter wave interactions 
and time-dependent noise in an essentially one-dimensional 



geometry. In contrast to previous work in the field of non- 
linear random waves |0 Holl . our initial condition does not 
correspond to a self-contained soliton because we assume re- 
pulsive interactions g > 0. Current experiments in wire traps 
have been hampered by the presence of a static field mod- 
ulation that leads to the fragmentation of the expanding atom 
cloud This is neglected here and makes a direct com- 

parison beyond the scope of our model. 

If we ignore spin flip processes for simplicity, magnetic 
noise in a microtrap above a planar substrate translates into 
a random potential with correlation function fl lTill 

(V{x,t)V(x',t')) « ~/6(t-t')C(\x-x'\) 



C(\x-x'\) = 



r- 



(2) 
(3) 



where 7 is the noise strength and the spatial correlation length 
i corr is of the order of the microtrap height iIeHi . If the poten- 
tial fluctuated only in time, 7 would correspond to the phase 
diffusion rate. Typically, 1/7 is a few seconds in /im sized 
traps 1 3, 4]. In the frequency range relevant for our model (up 
to fi r ), the noise spectrum is approximately flat 1 14]. 

For a single, typical realization of the noise, the evolution 
of the density n{x, t) = \ip(x,t)\ 2 according to Eq. is 
shown in Fig. \l\ A complicated fringe pattern appears due 
to the interference between the expanding condensate mode 
and the excitations generated by noise, with the fringe phase 
depending on the history of the noise. If we average over the 
evolutions in an ensemble of noise potentials, a smooth aver- 
age field ip c (x, t) = (ip(x,t)) with a decaying weight emerges 
(Fig. 0. This quantity would be revealed in a (as yet hypothet- 
ical) homodyne measurement of the Bose field, and we shall 
call it the 'coherent field' in the following. Note the analogy 
to the condensate order parameter in the symmetry breaking 
approach to Bose condensation when N c ~ N 

Another condensate definition, which is also applicable 
in U(l) covariant theories, is based on long-range order in 
the single-particle density matrix II 1611 . This quantity corre- 
sponds in our problem to the coherence function p(x, x',t) = 
(ip*(x, t)ip(x', t)). Its decay with increasing separation s — 
x' — x defines the coherence length l co h(t), a key concept of 
decoherence theory ll5i [l7ll . We introduce the spatial average 



(4) 



T(s,t) — / dx p(x, x + s,t), 
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FIG. 1: Expansion of a self-interacting Schrodinger field in a noisy 
potential, single realization. The normalized spatial density is plotted 
at different times given in the inset. Parameters in Eq.0: interaction 
strength gN — 10, noise strength 7=1, correlation length i corr = 
0.1. Harmonic oscillator units with respect to the initial confinement 
frequency are used: t i— » Sit, x i— > (Ti/mQ,) 1 ^ 2 x. The numerical 
solution uses a discrete space-time grid with time step dt = 0.1, and 
2 14 space points spaced dx = 0.0294 units. 




FIG. 2: Normalized density profiles of the coherent (noise-averaged) 
field for different expansion times. Symbols: Monte Carlo results, 
dashed lines: gaussian approximate solution to Eq. ^5}, \ t/j c (x, t) | 2 = 
{N/(y/2nu(t)) e _7t exp[-x 2 /{2u 2 (t)] with u(t) solving Eq. flOl . 
Inset: coherent fraction (relative particle number). Crosses (open 
circles): Monte Carlo results for 7 = 1 (7 = 0.1). Dashed (solid) 
line: exponential decay with decoherence rate 7 (renormalized rate 
7cff = 0.82 7). Units and all other parameters as in Fig.Q 

whose Fourier transform with respect to s is the momentum 
distribution, averaged over many realizations. This leads to 
IcohSp ~ 1 where dp is the width in momentum. The reduc- 
tion of the coherence length ('decoherence') is borne out in 
the results plotted in Fig. [5] Long-range coherence is also vis- 
ible: a fraction of the bosonic wave field is coherent across the 
full cloud size. We shall see that this fraction can be identified 
with the coherent field (ip(x,t)), reinforcing the analogy be- 
tween the condensate order parameter and the noise-averaged 
nonlinear Schrodinger field. 

In the two-component model of Bose-Einstein conden- 
sation, the condensate evolves according to a nonlinear 
Schrodinger equation including loss terms and interactions 
with the non-condensed density. The non-condensed com- 
ponent is described by a suitable kinetic theory Irj [Till . We 
have adapted this model to our problem by replacing the aver- 
age with respect to the field's density operator by the average 
over the evolutions in an ensemble of random potentials. For 
a similar approach, see fiH : the noninteracting case has been 
treated in 12011 . An essentially analytical solution has been 
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FIG. 3: Spatially averaged coherence function [Eq. 0] for differ- 
ent expansion times (values and symbol coding as in Fig. Sym- 
bols: Monte Carlo results, solid lines: kinetic theory [Eq. {8j] with 
Ti(s, 0) = 0, renormalized decoherence rate 7 e g = 0.82 7 and noise 
correlation length l e g = 1.25 i CO rr- Other parameters as Fig.Q 

obtained with two approximations: (i) we describe the fluctu- 
ations around the coherent field by the free space dispersion 
relation; (ii) we neglect in the nonlinear Schrodinger equation 
the interaction between the coherent field and its fluctuations. 
The resulting equations are 

id t A = - l -d 2 x i> c +g\Mx,t)\ 2 A ~ (5) 
for the coherent field, and the Boltzmann-type equation 

(d t +pd x )W i (x,p,t)= (6) 
J dp' S v {p - p') (W c (x,p', t) + Wi&p', t) - Wi(x, p, t)) 

for the Wigner representation W\ of the 'incoherent' field. 
This is the Wigner transform of the fluctuating part of 
the coherence function p\(x,x',t) = (tp* (x, t)tp(x', t)) — 
ip* (x, t)i/j c (x', t). In Eq. W c is the Wigner transform 
of il>*(x, t)il>r(x' , t), and the 'collision integral' involves the 
'cross section' |21] 

S V {P-P')=1 J d{X ~ n X,) C(x-x')e-^-^ x - x '\ (7) 

where C(x — x') is the normalized noise correlation func- 
tion Q. When Eq. is approximated by a Fokker-Planck 
equation and W c = 0, we essentially recover the decoherence 
model discussed by W. H. Zurek o]. 

The analytical solution to Eqs. (|5)|6l is derived using previ- 
ous results for a noninteracting gas 1 20, 22]. The new ingredi- 
ent is the nonzero average of the field that enters the collision 
integral as a source term. The basic idea of the analytical solu- 
tion is to perform a Fourier transformation with respect to both 
x and p, and to solve the resulting equation with the method 
of characteristics I23II . For the coherence function (|4j, this 
yields 

r(s,t) = r c ( s ,t) + r i ( s ,o)e-^[ 1 - c ^ 

t 

+jC(s) Jdrc-^-c^Tcis^-T). (8) 
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The coherent field contribution, T c (s,t), can be found approx- 
imately using a time-dependent Thomas-Fermi profile 1 24] or 
a gaussian ansatz 1I25I1 . We follow the latter method because it 
simplifies the Wigner transform and get 



r c ( S ,i)«JVe-^exp j-y 
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where the spatial width u(t) is the solution of 



+ u'(t))}, (9) 
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du \ 8 u 
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The effective interaction strength is g(t) = giVe -7 * / '(Ay/ir), 
and the initial condition minimizes U e g(u) = l/(8u 2 ) + 
g(t)/u + ^D, 2 u 2 where Q is the initial trap frequency. 

The result of the kinetic theory splits into a contribution 
with long-range coherence (the first term T c extends across 
the entire cloud size) and an incoherent part (the second and 
third terms). For the initial conditions considered here, the 
second term vanishes so that the incoherent fraction is not co- 
herent beyond the correlation length of the noise potential, as 
intuitively expected. (Recall that the spatial correlation func- 
tion C(s) [Eq. 0] decays on the scale £ C orr-) 

The simple kinetic theory outlined above captures quali- 
tatively the features observed in Monte Carlo simulations of 
Eq. Q, as shown by Figs.|2][3] We attribute deviations to the 
approximate treatment of interactions in the theory. 

The density profile of the coherent field is not exactly gaus- 
sian (Fig. |3 because for the chosen parameters, one already 
approaches the Thomas-Fermi regime. The coherent frac- 
tion of particles N c (t)/N = (l/N) Jdx \ip c (x, t)\ 2 , however, 
shows an exponential decay as predicted by Eq. (|3J. We find 
quantitative agreement with the kinetic theory when a 'renor- 
malized' scattering rate 7 c ff < 7 is used. Further simulation 
runs show that the ratio 7 e s /7 does not change significantly 
when the noise strength is reduced by an order of magnitude. 
A similar renormalization of the coherent field's loss rate has 
been observed in simulations of nonlinear pulse propagation 
with a randomly fluctuating phase mismatch 1 26] . An analyt- 
ical approximation for y e g /j is derived below. 

Similarly, the renormalization of the noise correlation 
length to ^ c ff > ^ corr reproduces quantitatively both the short- 
range and long-range behavior of the coherence function 
(Fig. [3}. The oscillations originate from the Thomas-Fermi 
density kink at the condensate border that is not captured by 
our gaussian ansatz. 

In a homogeneous condensate, atomic interactions suppress 
the generation of long-wavelength excitations by a spatially 
modulated potential. This is described by the vanishing of 
the dynamic structure factor S[k; n c ] in the limit k — > 1 15]. 
Let us consider for our problem the condensate as locally ho- 
mogeneous with a density n c (x) — |-0 C (^)| 2 - Including the 
structure factor in the scattering cross section, we suggest an 
improved approximation to the collision integral in Eq. (|6) 



no interactions 
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FIG. 4: Renormalized decoherence rate 7 e ff [Eq. <12H for in- 
creasing interaction strength. On the :r-axis is plotted the ratio 
icorr/Co oc ^/g, where £0 is the healing length at the trap cen- 
ter, £o = l/ \J 4:gn c (x — 0). Large dots: Monte Carlo data (see 
text). Small lower (upper) dots: numerical integration of Eq. 1 1 21 for 
a Thomas-Fermi (gaussian) density profile. Dashed lines: asymp- 
totics for a Thomas-Fermi profile at weak and strong interactions 
[Eq. 1141 1. thick solid line: interpolation 1151 . The decoherence rate 
is normalized to its value 7 in an ideal gas. 



S v (k)S[k; n c {x)]n c (x, t), 



(11) 



where Sv(k) [Eq. Q] is the wavevector spectrum of the 
noise, and the well-known Bogoliubov structure factor is lfl5ll 
S[k;n c ] = \k\/(k 2 + 4gn c ) 1 / 2 . In Eq. O we have re- 
placed the coherent field's momentum distribution by a <5- 
function, assuming a width u(t) ^> l corr . For the ideal gas, 
S[k; n c (x)] = 1, and we recover Eq. (|6}. In the general case, 
we find 



7eff(*) = 



N c (t) 



dx dk Sy (k)S[k; n c {x)]n c (x, t). (12) 



The relevance of interactions is now determined by the com- 
petition between the width l/l COII of Sy(k) and the width 
l/£(x) of the structure factor, involving the local healing 
length = (Agn^x))^ 1 / 2 . 

We have computed the integral fl!2i numerically for gaus- 
sian and Thomas-Fermi density profiles. Very similar results 
are found (see Fig.|4j when the reduction factor 7 e g /j is plot- 
ted versus the ratio Z CO rr/£o where the healing length £n is 
taken at the trap center. It turns out that the time-dependence 
only appears via Z C orr/£o oc 1/ yju[t). For the Thomas-Fermi 
profile, the x integral can be evaluated analytically, leading to 



7cff 3/ c 



7 



(13) 



where f(z) = z [z + (1 — z 2 ) arctanz] . Asymptotic analy- 
sis leads to the limiting behavior (dashed lines in Fig. |4jl 



1_0.88%^ forJ CO rr«£o, 

o c (14) 
7 1 3tt £ 
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for Co < 



d t Wi(x,k,t) = dk' S v (k')S[k';n c (x)]W c {x,k- 



t The ideal gas result y e g = 7 is recovered for £0 = oo> but the 
^ ' next order correction already comes into play for a moderate 
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interaction parameter g, since £o °c g 1 ^ 2 - Very strong inter- 
actions, where £o - > 0, significantly slow down decoherence 
within our model. In the intermediate regime £ ~ 'com the 
asymptotics dl4> is quite inaccurate, but the interpolation 



7 A + 8Z corr /(3<o)' 

with A w 1.15 provides good agreement. The decoherence 
rate extracted from the Monte Carlo results is also fairly well 
described by Eqs. Jl 3II15I . The two data points we have plot- 
ted correspond to two values for the healing length: based 
either on the numerically computed coherent field density or 
on its Thomas-Fermi approximation. 

Let us finally note that for large interactions or long ex- 
pansion times, the coupling to the incoherent fraction will no 
longer be negligible compared to the random potential and our 
approximate solution will lose accuracy. In this limit, we may 
expect that the noise-induced decoherence is replaced by an 
'open system dynamics' similar to the models discussed for a 
condensate at finite temperature 1 1 81 12711 . 

To summarize, we have evaluated the impact of weak mag- 
netic field fluctuations on coherent matter wave dynamics in 
atom chips. A quantum kinetic theory for a degenerate trapped 
boson gas subject to noise has been worked out in the mean 
field approximation and solved analytically, neglecting the 
backaction of excitations onto the coherent field. The compar- 
ison to numerical simulations demonstrates that interatomic 
interactions reduce the decoherence rate relative to the ideal 
gas. We have suggested an explanation in terms of the struc- 
ture factor of a quasi-homogeneous system that leads to an ac- 
curate agreement with the numerical data. Further investiga- 
tions will address the renormalization of the noise correlation 
length due to interactions and the impact of finite temperature 
in the initial conditions. 
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